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Hydrostatic gas distributions: global estimates of temperature and 
abundance 
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ABSTRACT 

Estimating the temperature and metal abundance of the intracluster and the intragroup 
media is crucial to determine their global metal content and to determine fundamental cos- 
mological parameters. When a spatially resolved temperature or abundance profile cannot be 
recovered from observations (e.g., for distant objects), or deprojection is difficult (e.g., due 
to a significant non-spherical shape), only global average temperature and abundance are de- 
rived. After introducing a general technique to build hydrostatic gaseous distributions of pre- 
scribed density profile in potential wells of any shape, we compute the global mass weighted 
and emission weighted temperature and abundance for a large set of barotropic equilibria and 
an observationally motivated abundance gradient. We also compute the spectroscopic-like 
temperature that is recovered from a single temperature fit of observed spectra. The derived 
emission weighted abundance and temperatures are higher by 50% to 100% than the corre- 
sponding mass weighted quantities, with overestimates that increase with the gas mean tem- 
perature. Spectroscopic temperatures are intermediate between mass and luminosity weighted 
temperatures. Dark matter flattening does not lead to significant differences in the values of the 
average temperatures or abundances with respect to the corresponding spherical case (except 
for extreme cases). 
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1 INTRODUCTION 

The amount of metals in the Intracluster Medium (ICM) and in 
the Intragroup Medium (IGM) gives us important clues about the 
past star formation activity of the stellar population of these galaxy 
systems, being it directly linked to the total number of supemovae 
exploded in the past and to the initial stellar mass function of the 
star formation epoch (e.g., Renzini et al. 1993). The metal content 
can also enlight how the enrichment proceeded, e.g., via stripping 
or galactic winds driven by SNe or AGN feedback, and has im- 
plications for both the ICM/IGM and galaxy evolution (e.g., Wu, 
Fabian, & Nulsen 2000; Finoguenov et al. 2001; Kapferer et al. 
2007). For these reasons the observational study of the metal con- 
tent of the ICM/IGM is growing fastly. After the first large com- 
pilation of (emission weighted) average abundance values of iron 
from EXOSAT, Einstein and GINGA observations (Arnaud 
et al. 1992), ASGA made metal measurements for many clusters 
(Fukazawa et al. 1994, Finoguenov et al. 2000, Baumgartner et al. 
2005). The average iron abundance was estimated to be 0.38±0.07 
and 0.21 ± 0.05 respectively for the cooling flow and non cooling 
flow clusters (Allen & Fabian 1998). In more recent times, the supe- 
rior quality of the XMM — Newton and Ghandra instrumenta- 
tion has allowed for more accurate determinations of the elemental 
abundance pattern (e.g., Tamura et al. 2004, Fukazawa et al. 2004, 



Durret et al. 2005, Sanders & Fabian 2006, de Plaa et al. 2007, 
Finoguenov et al. 2007, Rasmussen & Ponman 2007). Nowadays, 
these studies are carried on also with Suzaku (e.g., Matsushita et 
al. 2007, Sato et al. 2007). 

Similarly to the metal abundance, the hot ICM/IGM tempera- 
ture is also one of the most important and commonly used global 
observables: it is used as a proxy for the total mass of the system 
(e.g., Voit 2005), from which the clusters can be used as probes for 
fundamental cosmological parameters (e.g., Henry & Amaud 1991, 
Henry 1997, Nevalainen et al. 2000, Amaud et al. 2005). Temper- 
ature profiles have been built with improved quality in the recent 
past (e.g., Arnaud et al. 2005, Pointecouteau et al. 2005, Vikhlinin 
et al. 2005, 2006, Pratt et al. 2007, Rasmussen & Ponman 2007). 
Since the ICM/IGM are not isothermal, ideally the mass weighted 
temperature should enter the computation of quantities to be used 
for cosmological tests. 

From a more quantitative point of view, the amount of the 
mass of metals in the ICM/IGM is given by 



= / p(x)Z(x)d'x, 



(1) 



where p and Z are the true three dimensional gas density and abun- 
dance profiles. Thus, the mass weighted average abundance is given 
by 
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<z>= 



(2) 



where Afgas ~ J p{x)d^x is the total hot gas mass. Similarly, the 
mass weighted average temperature is 

/ p(x)r(x)d3x 



<T> = 



(3) 



Unfortunately, there are at least three serious problems with esti- 
mating <Z> and <T> from observations: 1) for many clus- 
ters/groups we do not know the intrinsic shape of the gas distri- 
bution and the viewing angles under which we are observing it; 
therefore, one cannot uniquely deproject observed quantities (ob- 
tained in general from X-ray data) to derive p, T, and Z; 2) even 
for spherically symmetric systems, deprojection is a demanding nu- 
merical process, very sensitive to the properties of the instrumen- 
tal PSF and to measurement errors (e.g., Finoguenov & Ponman 
1999); 3) in many cases only a single spectrum can be extracted 
for the whole gas, and only an average abundance and temperature 
can be obtained; this happens when there are not enough counts for 
a spatially resolved spectroscopy, e.g., for distant clusters/groups 
(Hashimoto et al. 2004, Maughan et al. 2007, Baldi et al. 2007 for 
recent observations with Chandra and XMM-Newton). In partic- 
ular, the average abundance and temperature mentioned in point 3) 
above are not those given in eqs. ([2l(-([3l(, but are in practice lumi- 
nosity weighted quantities (e.g., Mathiesen & Evrard 2001, Maz- 
zotta et al. 2004, Maughan et al. 2007, Rasia et al. 2005, Kapferer 
et al. 2007) that can be defined as 



<Z>L = 

and 

<T>i= 



/ Ex(Ci,e2)V(Ci,6)rf6rfC2 



/ Sx(gi,g2)rp.(gi,C2)dCirfC2 

Lx 



(4) 



(5) 



where (^1,^2) are the coordinates of the projection plane. Ex is 
the X-ray ICM surface brigthness, Zp,- and Tpr are the luminos- 
ity weighted projected abundance and temperature, and Lx = 
J I2xd£,id£^2 is the total X-ray luminosity (see Appendix Al). 

It is then natural to investigate the relation between the quan- 
tities in eqs. (2)-(3) and (4)-(5). For example, Rasia et al. (2008), 
using mock XMM — Newton spectra for a sample of simulated 
clusters, find that the iron abundance inferred from such spectra is 
very close to the projection of the emission weighted values of Z 
(i.e., <Z>l), at least for thermal components of kT > 3 keV and 
kT < 2 keV. Kapferer et al. (2007), again using simulations, simi- 
larly find that for kT > 3 keV the X-ray emission weighted abun- 
dance is close within few percents to that derived from the analy- 
sis of synthetic X-ray spectra. Unfortunately, neglecting a possible 
spatial variation of the metal abundance can lead to largely wrong 
estimates of when using <Z>l instead of <Z> in eq. {2} 
(Amaud et al. 1992). In fact, iron distributions peaked towards the 
cluster/group center have been revealed in many cases (Fukazawa 
et al. 2000, Ettori et al. 2002, Sanders & Fabian 2002, Matsushita 
et al. 2003, Bohringer et al. 2004, Tamura et al. 2004). Motivated 
by this, in an exploratory study Pellegrini & Ciotti (2002) showed 
that in these cases <Z> can be significantly smaller than <Z>l- 
Successively, De Grandi et al. (2004) confirmed this result for their 
sample of cooling core clusters, for which they estimated <Z> to 
be ~ 15% smaller than <Z>l- 

It is also well accepted that the ICM/IGM have a temperature 
structure that was established by gravitational and non-gravitational 



processes, as radiative cooling and heating by active galactic nu- 
clei (see Borgani et al. 2005, Vikhlinin et al. 2005, Piffaretti et al. 
2005, Amaud et al. 2005, Donahue et al. 2006). Efforts have been 
made recently to understand the meaning of the temperature de- 
rived from spectroscopic observations when the ICM/IGM has a 
complex thermal structure (Mazzotta et al. 2004, Rasia et al. 2005, 
Vikhlinin 2006, Nagai et al. 2007). Mazzotta et al. (2004) found 
that the observed temperature, recovered from a single tempera- 
ture fit to the spectrum of a plasma with components at different 
temperatures (but all continuum-dominated, i.e., with kT>3 keV) 
and extracted from Chandra or A/— Newton data, is well ap- 
proximated by a "spectroscopic-like temperature" T^i (see Sect.[3](. 
Vikhlinin (2006) extended this previous work and proposed an al- 
gorithm to accurately predict Tsi that would be derived for a plasma 
with components in a wider range of temperatures (kT>0.5 keV) 
and arbitrary abundances of heavy elements. From the analysis of 
mock spectra of simulated clusters, it was found that Tsi is lower 
than the emission weighted temperature <T>l, with consequences 
for using the observed M — T relation to infer the amplitude of the 
power spectrum of primordial fluctuations (Rasia et al. 2005). 

Here, extending the preliminary discussion of Pellegrini & 
Ciotti (2002) based on spherical models, we estimate how much 
discrepant <Z>l and <Z>, and <T>l (or Tsi) and <T> are, by 
using different plausible profiles for p, T and Z obtained assuming 
hydrostatic equilibrium within triaxial mass distributions resem- 
bling real systems. In particular the models are constructed by using 
a technique that allows for building analytical barotropic gas distri- 
butions with prescribed density profiles departing from spherical 
symmetry. These new models extend the class of equilibria usually 
considered in the literature beyond isothermal or polytropic models 
(i.e., Suto, Sasaki & Makino 1998; Pellegrini & Ciotti 2002; Lee & 
Suto 2003, 2004; Ostriker, Bode & Babul 2005; Ascasibar & Diego 
2007). In the computation of the averages, our approach takes also 
advantage of the Projection Theorem, from which it follows that 
<Z>L and <T>L are independent of the specific direction of 
the line-of-sight, and can be calculated using the intrinsic three- 
dimensional quantities of the models, with a much easier procedure 
that avoids projection and surface integration. 

The paper is organized as follows. In Section 2 we present the 
models of the dark matter halos and the procedure to build fully 
analytical hydrostatic configurations in potentials of triaxial shape, 
for gas distributions corresponding to truncated quasi-isothermal 
models, quasi-polytropic models and modified /3 models. In Sec- 
tion 3 we describe the results and in Section 4 we summarize the 
main conclusions; technical results are reported in the Appendix. 



2 THE MODELS 

2.1 Density profiles for the gravitating mass 

The density of the (dark) mass distribution is the generalization to 
the triaxial case of the so-called 7-models (Dehnen 1993, Tremaine 
et al. 1994): 



Af(3 - 7) 



^ 47rr-?(l-e)(l-r?) m^(l + m)4-T' 
where 



2 ^ , 



y 



+ 



(l-e)2r| (l-r7)2r, 



(6) 



(7) 



M is the total dark mass of the system, is a characteristic 
scale, and the pair (e, rj) parameterizes the flattening along the y 
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and z axes respectively. The mass distribution is spiierically sym- 
metric when e = ?7 = 0, and M remains constant for different 
choices of the flattening. For simplicity, we restrict to the 7 = and 
the 7 = 1 cases: in the former, the density profile shows a central 
"core", while in the latter the Hernquist (1990) profile is recovered 
in the spherical limit. Note that the 7 = 1 models have the same ra- 
dial trend, in the central regions, as the profile obtained from high 
resolution cosmological simulations (Dubinsky & Carlberg 1991; 
Navarro, Frenk & White 1996), while they are steeper at large radii 



(oc r 



instead of oc r ). Even though not required by the tech- 



nique described in Sect. 2.2, in our analysis we used the potential 
profiles obtained by means of homeoidal expansions of the true po- 
tential at fixed total mass (e.g., Muccione & Ciotti 2003, 2004; Lee 
& Suto 2003, 2004; Ciotti & Bertin 2005, hereafter CB05). This 
approach has the advantage of avoiding the numerical integration 
needed to recover the potential (e.g. Binney & Tremaine 2008), and 
the formulae obtained are a very good approximation of the exact 
potential associated with eq. l[6). 

Homeoidal expansion applied to the 7 = model shows that 



2r + 1 



2(r - 
121n(r 



3r^ + 12r- + 8 41n(r + l) 
2(r + 1)2 r 



e + 77 



3r^ + 22r^ + 30r + 12 



2 I 2 
ey + jyz 



(r + l)3 

where 4> = GM4>/rc, and the value of the central potential is (jio 
— (3 + e -|- 77) /6. For the 7 = 1 model 

\ + 2 21n(r + l)" 



(8) 



r + 1 



61n(r+l) 2r2 + 9r + 6 



(r + iy 



2 1 2 

ey + rjz 



(9) 



and (po = —(3 + e + rj)/3. In the formulae above the radial co- 
ordinates are normalized to Vc, and in both cases the expansion 
holds for 1 > Srj — e (see Appendix A in CB05). Thus, in prin- 
ciple the maximum deviation from spherical symmetry is obtained 
for r] = e — 0.5, corresponding to a prolate system of axis ra- 
tio 2:]13]. Finally, the virial temperature of the system (defined as 
SkMTvir = \U\, where U is the gravitational energy) in the limit 
of small flattenings, and independently of the specific density pro- 
file g(m), is given by 



GM firriu 



(10) 



where ^ is the mean particle weight, mu is the proton mass, k is the 
Boltzmann constant and r-vir is the virial radius of g in the spherical 
limit (Muccione & Ciotti 2004). Here r^ir = lOr-c and 6rc for the 
7 = and 7 = 1 models, respectively. Note that, for fixed M and 
''vir, Tvir lucreases for an increasing flattening. 

Summarizing, the potential is determined by assigning the two 
flattenings e and rj, and by choosing the mass M, the slope 7, and 
re- The latter step is done via the relation rvir ~ rvir(M) holding 
for dark matter halos obtained from cosmological simulations in a 
flat ACDM cosmological model {Qm = 0-3, = 0.7, h = 0.7, 
where the Hubble constant is defined as lOO/i km s^^ Mpc^^), 
as derived, e.g., by Lanzoni et al. (2004). For example, for a mass 
M = 3.5 X lO^'^h'^MQ we adopt rvir = 1.4/i"^ Mpc, so that 



^ As shown in CB05, for large flattenings the expanded density deviates 
from an ellipsoid, being more similar to a toroid; however the shape of the 
equipotential surfaces is very similar to that of ellipsoidal systems. 



Tc = O.Uh^^ Mpc for the 7 = model, and = 0.23ft~^ Mpc 
for the Hernquist model, with Tvir ~ 2.3 keV (spherical case). For 
M = 1.0 X 10^^/i"^A/q, r-vir = l.8h~^ Mpc and Tvir = 5.1 keV 
(spherical case). We also derived the commonly used r2oo and rsoo 
radii (within which the average mass density is respectively 200 
and 500 times the critical density at redshift zero for a flat ACDM 
cosmological model). Independently of 7 = or 7 = 1, r2oo — 
0.7r-vir and rsoo ~ 0.5rvir for M = 3.5 x lO^'^ft-^Af©, and 
r2oo ==; O.Srvir and rsoo 0.6rvir for M = 1.0 x 10^^/i"^Mq. 
Remarkably, the ratios r2oo /rvir and rsoo / '"vir are very similar to 
those typical of the Navarro et al. (2006) profile of same total mass 
and virial radius. 



2.2 The hydrostatic equilibrium models 

Once a dark matter distribution is chosen, we build hydrostatic 
equilibrium models for the gas within it, assuming that the gas 
mass does not contribute to the gravitational field, and that the gas 
is perfect so that its pressure is p — kpT/ ^mn- Our procedure is 
based on the well known result that pressure, density and temper- 
ature in hydrostatic equilibrium are all stratified over isopotential 
surfaces (e.g., Tassoul 1980]3- In other words, hydrostatic config- 
urations are barotropic, i.e. p = p{p), which allows us to solve 
the hydrostatic equation Vp = ~p\7(f> for potentials of general 
shape. Therefore, the method is fully general: the only additional 
simplifying assumption is that the potential has a finite minimum 
00 at the center and vanishes at infinity. With this method we could 
also study the effect of substructures by superimposing different, 
off-centered dark-matter halos. 



2.2.1 Truncated quasi-isothermal models 

The following is a family of exact equilibria that generalizes the 
classical isothermal models 



P = Po e 



(j}—<j}Q 



kTo 



(11) 



where p is the isothermal equilibrium stratification of temperature 
To in a generic potential 0, and (po and po are (for example) the 
central potential and the central gas density. As usual for isothermal 
equilibria the total mass diverges, and a truncation surface (outside 
which p = 0) must be introduced. This should be done preserving 
the barotropicity of the distribution. In practice, the truncation sur- 
face must be an isopotential surfac^f]- In addition, to avoid unphysi- 
cal density jumps, it is natural to truncate the system by subtracting 
to eq. i ll It (the parent distribution), its value on some isopotential 
surface (j)t, so we consider the new density distribution 



p = Po e "0 



<l>o / _ <t> 
■3o (e ^ 



■e W 



(j) < (t>t, 



(12) 



while the quasi-isothermal equilibrium temperature associated with 
eq. ( 112b is obtained from eq. (A6) as 



If ^ vai'ies, it is actually the ratio T//i to be stratified over the isopotential 
surfaces, but here we neglect the very small p variations due to the adopted 
abundance gradients. 

Note the analogy with stationary truncated stellar systems where, accord- 
ing to the Jeans theorem, the truncation surface must be defined in terms of 
the isolating integrals of the motion. At the truncation surface, the normal 
component of the velocity dispersion tensor (the temperature analogous) 
vanishes (e.g. Ciotti 2000). 
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Figure 1. Radial profiles of density and temperature for truncated quasi-isothermal spherical gas models in the Hernquist potential (left, rt = rvir = 6rc) and 
in the 7 = potential (right, rt = r^\^ = lOrc), with T^ir = 2.3 keV. Three values of Tq of the parent isothermal model have been chosen: Tq = 0.8 T^ir 
(solid lines). To = 1.5Tvir (dotted lines) and To = STvir (dashed lines). The distance from the center r is normalized to r200 as defined at the end of 
Sect.[2Jl 



T_ 

To 



= 1 



<t>t-4 

e '^0 



- 1 



(13) 



A different approach, that we do not explore here (but that 
could be easily implemented in our scheme), would be that of fix- 
ing the pressure to some prescribed value on the truncation surface, 
by imposing a finite density jump at 0t, as done in Ostriker et al. 
(2005). Note that the central values of T and p of the truncated 
distribution are not po and To of the isothermal parent distribution 
in eq. i ll U . and the temperature at the truncation surface vanishes. 
Formally, the untruncated case (i.e., the true isothermal case) is re- 
covered for (/!>t 00, or for To —> 0. At the opposite case, i.e., for 
very large To, the following asymptotic behavior is obtained: 



p ~ po- 



l3o > 1. 



(14) 



f3o ' 2/3o 
In this limit the temperature distribution becomes independent of 
To, and p cc p^. Also the asymptotic density profile, for an assigned 
gas mass, is independent of Tq. 

Summarizing, a quasi-isothermal model is determined by 
choosing a mass model as described in Sect. 2.1, and by assum- 
ing (pt ~ (t>{rvii) (that we arbitrarily fix along the a;-axis, see 
eqs. |[8)-(9)). Then a To is chosen and po is obtained by imposing 
that the total Mgas of the truncated distribution equals a prescribed 
value. Figure[T]shows the density and temperature profiles of quasi- 
isothermal equilibria in a 7 = 1 and 7 = spherical mass distri- 
bution. The total dark matter mass is M = 5 x lO'^ Mq and we 
assume Mgas = 0.14M, according with the direct measurements 
of gas mass fractions of LaRoque et al. (2006), for the concordance 
flat ACDM model. As expected, flatter temperature profiles are ob- 
tained for lower values of To /Tvir, while for high values of To /Tvir 
the density profile tends to the limit distribution l ll4t . In case of in- 
termediate dark matter flattenings (e.g., e = 0.1, 77 = 0.3), the 
maximum flattening of the gas distribution is ~ 0.10 in the {x, z) 



plane, while in the e — rj = 0.5 case the maximum gas flattening is 
~ 0.16. These figures are similar in the 7 = and 7=1 models, 
and go in the expected direction. The reason for this lies in the well 
known fact that the gas density and temperature distributions are 
stratified on equipotential surfaces, that are much less flattened than 
the mass distribution that produces them (e.g., Binney & Tremaine 
2008). Therefore, even for the flattest mass distributions that can be 
allowed, the corresponding density profiles keep roundish. 

2.2.2 Truncated quasi-polytropic models 

Polytropic models are equilibrium statifications for which p = 
Po{p/po)^ and T = To{p/ po)^^^ , with the polytropic index 
1 < r < 5/3, and po and To are (for example) the central val- 
ues of the gas density and temperature, respectively. These models 
are more complicate than isothermal stratifications. In fact, in this 
case the solution of the hydrostatic equilibrium can be written as 

l3o= (15) 



_ T _ r - 1 

^ rpo 



where /3o now refers to the central value of the temperature. It fol- 
lows that, given the depth of the potential well, a critical tempera- 
ture 



T 

^ Cl 



r — 1 pmH|<^ol 



(16) 



r k 

exists so that for To > Tcr the distribution in eq. l ll5t is untrun- 
cated, and the total gas mass diverges. For To = Tot < Tcr instead 
a truncation value (jjt defined by the identity 

r - 1 pmH{(f)t - 0o) 
= — k ^^^^ 

exists, so that r((?it) = 0. Alternatively, having fixed the two val- 
ues > ^t > 00 for the potential, only one temperature Tot 
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Figure 2. Radial profiles of density and temperature for truncated quasi-polytropic spherical gas models of index F = 1.2, in the Hemquist potential (left) 
and 7 = potential (right), for the same T^\^ = 2.3 keV and mass parameters of Fig. [T] The central temperature of the parent polytropic distribution is 
To = 2.6Tvir (solid lines) and To = 4.0Tvir (dotted lines). 



exists that produces a naturally truncated polytrope at the surface 
(j) = (j)t. However, it can be useful to have a whole family of quasi- 
polytropic models truncated at (j)t for all temperatures To > Tot . 
This can be obtained following the same approach as in Sect. 2.2.1. 
Thus, for given <jit and To > Tot, we introduce the truncated den- 
sity 



(18) 



po \io/ \io, 

where T is the temperature of the parent model l llSb . and Tt is 
its value at (j>t\ of course Tt = for To = Tot. Following the 
method described in Appendix A, the quasi-polytropic equilibrium 
temperature corresponding to eq. dlSt is 



T_ ^ \To) VTo 
To 



3o \To) 



Tt\ — 



(19) 



where the temperature distribution at the r.h.s. is that given by 
eq. J15l >. Summarizing, after having choosen a dark matter distribu- 
tion and the value (j>t = (/"(^vir) as in the quasi-isothermal case, the 
associated Tot is calculated. A truncated quasi-polytropic model is 
then determined by fixing a temperature To > Tot, so that Tt is 
determined through eq. (15), and po is obtained so that A/gas of the 
truncated distribution (18) coincides with the required value. 

We remark that the pair ([T8ll-([T9ll when To = Tn reduces to 
the polytrope naturally truncated at <j)t, while for very high values 
of the central temperature 



p ~ po- 



rpo 



2Po ' 



A) > 1, 



(20) 



and, as in the quasi-isothermal case, the temperature distribution 
becomes independent of To. For reference, from eqs. MA\ and 
(20) it follows that for the limit 7 = 1 models the ratio of the true 



central gas temperature r(0) to Tvir is ~ 7.7, while in the limit 
7 = models it is T(0) ~ 6.2rvir. 

Figure[2] shows the density and temperature profiles for quasi- 
polytropic spherical models with F = 1.2 (a value reported to pro- 
duce a good fit of some observed temperature profiles for the ICM, 
Markevitch et al. 1998) in the same potentials adopted for Fig.[T] As 
for the truncated quasi-isothermal models, steeper density profiles 
in the central regions are obtained for the 7 = 1 than for the 7 = 
potential, to balance the steeper potential well (even though in the 
quasi-polytropic case the steepening can be minor, being in part 
compensated by the temperature increase towards the center). Note 
that models analogous to the "coldest" quasi-isothermal models in 
Fig.[T]do not exist because from eq. Ml\ the minimum admissible 
temperature Tot is 2.irvir for 7 = and 2.&Tvii for 7 = 1. As in 
the quasi-isothermal cases, also here the effect of dark matter flat- 
tening on the density and temperature distributions is quite modest. 
In fact, being the gas stratified on the potential, the flattenings of the 
gas distributions are the same as described at the end of Sect. |2.2.T| 



2.2.3 Truncated modified (3 models 



The models introduced in the previous Sects. [ZZTI and 12.2.21 are 
just two special barotropic families built starting from prescribed 
relations p(p); as a consequence, their density profile is somewhat 
out of control. Here we show how to derive the temperature dis- 
tribution for an hydrostatic gas of assigned density profile in an 
external potential well deviating from spherical symmetry. We call 
this approach "density approach'0 and technical details are given 
in Appendix A2. In practice, the idea behind the method is to con- 
struct the spherical barotropic solution for a given gas density pro- 



* For the more complicate case of the construction of rotating, baroclinic 
gaseous distributions, see Barnabe et al. (2005). 
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Figure 3. Gas density and temperature radial profiles for a TMB model (with f} = 2/3) in equilibrium within the spherical Hernquist potential used in Figs. 
1 and 2. From left to right the panels refer to a = 0, a = 1 and q = 2. In each panel the solid, dotted and dashed lines refer to Tg/rc = 0.4, 1, and 1.6, 
respectively. Note that consistently with eq. I2U the profiles of the a = 2 models are independent of rg/rc. 



file in a given spherical potential, and then to deform (maintaining 
the equilibrium) the potential and the gas density distribution: this 
is accomplished by constructing the integral function H . 

As relevant case for the present discussion, the starting density 
distribution is a spherical truncated modified /3-model (hereafter 
TMB) 

for r < rt, with rg a core radius and rt a truncation radius. This 
density profile is a modification of the well known /3-model (Cav- 
aliere & Fusco Femiano 1976) and its generalization (Lewis et al. 
2003). In particular, the density is proportional to for r ^ 0, 
and to r~^^ for rg ^ r ^ rt. A finite gas mass is obtained for 
< a < 3, and for 13 > 1 no truncation would be required. 
Here the introduction of rt is needed because 0.5</3<0.8 from fits 
to observed ICM profiles (e.g., Mohr et al. 1999, Jones & Forman 
1999). For a spherical Hernquist potential, the density approach ap- 
plied to eq. l l21b leads to the function 

p(^) 

poh^P ^ (1 - + 62)vl-2 _ + l](3/3— )/2 

1 (2T) 

(l-*t)'^[(l + b2)4'2-2*t + l](3'5-°'/2' ^ > 

where '5 = (jy/i'o is the Hernquist potential normalized to its cen- 
tral value, and b = r^/rc- Note how the two limiting cases of very 
small and very large 6 correspond to truncated power law gas dis- 
tributions: p oc r"^'' — r'^^^ for 6 — > 0, and p oc r~" — r^° 
for b — > oo. The function H{'i/) needed to determine the tempera- 
ture distribution (eqs. [A4], [A6]) cannot be expressed in terms of 
elementary functions for generic values of a and /3; however, sim- 
ple cases are obtained for a = 0, 1, 2 and /9 = {a + n)/S with 
n non-negative integer. The explicit formulae for a = 0, 1, 2 and 
/3 = 2/3 (that falls within the observed range quoted above) are 



provided in Appendix A2, and hereafter only these values will be 
used. Thus small values of rg/rc correspond to models converging 
to the truncated r~^ profile, independently of the specific value of 
a, while for a = 2 the distribution is independent of b. The final 
step of the procedure is to substitute the deformed potential given 
in eq. ^ in eq. i22i and in the function H, since by construction all 
the resulting formulae are still exact when the potential is deformed 
to the axisymmetric or triaxial case. 

Figure [3] shows the density and temperature profiles for the 
a = 0, 1, 2 spherical cases. The temperature decline in the q = 1 
and a = 2 models compensates the steep increase of p, in order 
to produce the pressure gradient needed to balance the imposed 
gravitational field. In a broad sense, this behavior is similar to that 
of the velocity dispersion profile in the central regions of isotropic 
Hernquist or R^^"' models (Ciotti 1991). For a = 1, lower val- 
ues of b correspond to a more important central peak of the density 
profile and a more important decline of the temperature in the cen- 
tral region. Thus, although the central temperature drop is not due 
to cooling, these models provide an interesting phenomenological 
description of cool-core systems. The opposite behavior is shown 
by the a = Q models, in which the flat-core gas density requires 
central temperatures higher than in all the other cases. Finally, the 
introduction of flattening in the dark matter halos does not lead to 
significant deformations in the gas density distributions, with max- 
ium deviations as reported at the end of Sects. [Z2n and l2?2.2l 

2.2.4 Comparison with observed ICM properties 

Even though the aim of this work is not to construct models repro- 
ducing in detail the observed ICM properties (which is hard within 
the simple framework of hydrostatic equilibrium of single-phase 
gas in smooth potential wells), we briefly comment here on how 
the obtained equilibria compare with observations. In general, the 
quasi-isothermal and polytropic models, and the TMB models with 
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Figure 4. The mass weighted <T> (black), luminosity weighted <T>£ (red) and spectroscopic-like T^i (green) temperatures for the truncated quasi- 
isothermal models (left) and the truncated quasi-polytropic models with F = 1.2 (right). Solid lines refer to the 7 = potential, dashed lines to the Hernquist 
potential (see Sect. l3.1l for more details). 



a = 0, are similar to "non cool-core" systems, while the a = 1 and 
Q = 2 TMB models, where the temperature profile is decreasing 
towards the center, are similar to "cool-core" systems. 

In the family of non cool-core models, quasi-isothermal dis- 
tributions can be built with arbitrarily low temperatures, becoming 
more and more similar to the standard isothermal models. Quasi- 
polytropic models instead, once the truncation potential is fixed, 
cannot be built with a central temperature smaller than a limit tem- 
perature roughly corresponding to the depth of the dark matter po- 
tential well. In the past, polytropic models with T — 1.2 have 
been used to reproduce the external regions of ICM observed (e.g., 
Markevitch et al. 1998, Piffaretti et al. 2005) and simulated (Os- 
triker et al. 2005). 

In the family of cool-core models, TMB distributions with 
Q = 1 or a = 2 show temperature profiles in good agreement with 
those observed by Chandra and XMM — Newton (e.g., Allen 
et al. 2001, Kaastra et al. 2004, Vikhlinin et al. 2005, 2006); on av- 
erage, these profiles reach a maximum near and then decline at 
larger radii, reaching ~ 0.5 of their peak value near r ~ 0.5rvir. 
In addition, not only the profile shapes of these TMB models 
are similar to the observed ones, but also their temperature val- 
ues, when rescaled to the mass-weighted temperature within rsoo 
(Tsoo), agree with the observed values (as those shown by Viklinin 
et al. 2006). The relation between <T> and Tsoo will be briefly 
addressed at the end of Sect. 3.1. 



2.3 The abundance profile and the emissivity 

In addition to the dark matter potential well and the hydrostatic gas 
distribution, the third ingredient of our models is the metal distribu- 



tion. In the numerical code the metal abundance profile is assumed 
to be stratified according to a formula which generalizes to the el- 
lipsoidal case the observed abundance profiles (Ikebe et al. 2003, 
De GrandElet al. 2004, Vikhlinin et al. 2005), i.e. 

7 = ^0 2 _ £^ 2^ 

where the central metallicity is Zq = O.SZq, the slope = 0.18 
and the metallicity scale-length rz = 0.04r-vir. In addition, the flat- 
tening of the metallicity distribution is the same used for the dark 
matter distribution. Obviously, we are not attaching any special 
physical reason to this last assumption, except to have flatter metal 
distributions in flatter systems, and to reduce the parameter space 
dimensionality. In any case, we also explored cases where the met- 
als are stratified exactly on isodensity surfaces, i.e. Z oc (p/po)^, 
without finding significant differences with the case of eq. | |23] |. 
The emissivity adopted in the code is given by 

£ = n^nnAiT, Z) (24) 

where and nu are the number densities of electrons and hydro- 
gen. The cooling function A(r, Z) has been calculated over the en- 
ergy interval 0.3-8 keV with the radiative emission code APEC for 
hot plasmas at the collisional ionization equilibrium (Smith et al. 
2001), as available in the XSPEC package (version 12.2.0) for the 
solar abundance ratios of Grevesse & Sauval (1998). With APEC 
we have computed a matrix of values for A(r, Z) for a very large 

^ The De Grandi et al. (2004)'s formula does not have the square of the 
radial coordinate, but we found that the square is needed to match the data- 
points in their Fig. 2. 
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Figure 5. The mass weighted <T> (black), luminosity weighted <T>i 
(red) and spectroscopic-like T^i (eq. 1291 . green) temperatures calculated 
for TMB gas models within the 7 = 1 potential, for a = (dotted lines), 
o = 1 (dashed lines) and a = 2 (solid lines), as a function of r^/rc- 
For a = 2 the average temperatures are independent of Tg/rc, consistently 
with eq. (21) with /3 = 2/3. For a = the value of <T>l /T^ir is 
displayed only up to rg /rc = 1 since for higher rg / rc the central temper- 
atures become very high (Fig.|3). 



set of temperatures and metallicities. Note that the cooling function 
can be written as 



A(r, Z) =A(r,0) [l + Zg{T,Z)], 



(25) 



where Z is in solar units, A(r, 0) is the function in the case of no 
metals, and 



9{T,Z) 



A(r, Z) - A(T, 0) 
ZA{T,0) 



(26) 



It turned out that the function g is almost exactly independent of 
Z, so that eq. i25\ with g = g{T) exploits the nearly perfect lin- 
ear dependence of the function A(r, Z) on abundance. In order to 
speed up the numerical code we computed non-linear fits of the 
functions A(r, 0) and g{T) valid over the temperature range 0.1- 
16 keV (with maximum deviations from the APEC values < 1%) 
and reported in Appendix B. We remark that g{T) declines steadily 
with increasing T, from ~ 42 down to ~ 0.14, which has the con- 
sequence that for high values of Zg{T) then A(T, Z) oc Zq and 
Lx oc Zo, while for low values of Zg{T) both A(T, Z) and Lx 
are independent of Zo- 



3 RESULTS 

For each model the quantities <Z>l and <T>l are not com- 
puted through projection, but directly as volume integrals. In fact, 
from the Projection Theorem (Appendix Al) eqs. (4)-(5) can be 
also written as 



<Z>L = 

and 

<T>L = 



/ £:[p(x),r(x),z(x)]Z(x) d^x 



Lx 



/ f[p(x),r(x),Z(x)]T(x)d3x 



(27) 



(28) 



where £ is the emissivity in the 0.3-8 keV band due to gas in the 
temperature range 0.1-16 keV. As anticipated in the Introduction, 
for each model we also compute the spectroscopic-like tempera- 
ture. Following Vikhlinin (2006), this is estimated as 



Tsi = arTcont + (1 — a;)T'iine 



(29) 



where Tcont and Tunc are the continuum-based and line-based tem- 
peratures for the composite spectrum and a; is a parameter that mea- 
sures the relative contribution of the line and continuum emission 
to the total flux. To evaluate Tcont, riinc and x, three functions of 
the temperature are needed; these depend on the instrument in use, 
energy band, redshift and neutral hydrogen absorbing column. As 
a representative case we chose to simulate observations made with 
the Chandra CCDs over the 0.3-8 keV band, for a plasma at zero 
redshift and zero absorbing column. A. Vikhlinin kindly provided 
us with the tabulated values of the required functions, that we fit- 
ted with the same high precision method described in Appendix 
B and we then inserted in our code. We recall that the method of 
Vikhlinin (2006) holds for thermal components of fcT>0.5 keV. 
We also computed the estimate of Tsi proposed by Mazzotta et al. 
(2004) for plasma components at fcT>3 keV: 

/p2(x)r^-l/2(x)d3x 



T,i(5) 



/ p2(x)T'5-3/2(x)d3j 



(30) 



where S = 3/4 for observations obtained with Chandra and 
XMM - Newton. 

It is important to note that if Z and T do not depend on x, 
thei{f| <Z>L=<Z>= Z and <r>L=<r>= T^i = T, even 
for density distributions depending on x; otherwise the variously 
weighted quantities differ, in a way dependent on the spatial distri- 
bution of p, T and Z. A quantitative estimate of these differences 
is the task of the following Sects. [3T| and |3. 21 

For each model all integrals have been calculated numerically 
with a double-precision code. The integration scheme employs a 
linear interpolation of the gridpoint-defined variables, and the num- 
ber of grid points in the positive octant of the [x ,y ,z) space is 
X Tij, X Tiz = 300 X 300 X 300. Checks of the code have been 
performed by calculating (with both linearly spaced and logarith- 
mic grids) the total masses of strongly peaked triaxial distributions 
whose values are known analytically, and also mean value temper- 
atures of special distributions for which the expected values can be 
calculated analytically (see Section 3.1), obtaining errors < 0.1%. 



3.1 Temperature averages 

Figures |4] and |5] show the trend of the mass weighted tempera- 
ture <r>, of the luminosity weighted temperature <T>l and 
of Tsi (eq. 1291 ), as a function of To /T^ir in the quasi-isothermal 
and quasi-polytropic cases, and of r^/rc in the TMB case. As for 
Figs. lU3l the gravitating mass is a spherical 7 = or 7 = 1 model 
with M = 5 X 10^"' Mq - also the range of To/Tvir and r^/r^ is 
the same used for Figs. lU3l 



^ This identity in the case of T^i of Vikhlinin (2006) can be easily proven 
considering the explicit form of the functions entering in eq. j29t . 
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Figure 6. The mass weighted <Z> (black) and luminosity weighted <Z>]^ (red), calculated for quasi-isothermal equilibria (left) and for quasi-polytropic 
equilibria with F = 1.2 (right), within the 7 = (solid lines) and 7 = 1 (dashed lines) potentials. The abundance profile is that of eq. (23), with Za = O.SZq. 



A first general result is that at this mass M the two T^i esti- 
mates of eqs. l |29l l and l l3Qt agree within ~ 10% for all the explored 
models. The reasons for this are the relatively flat shape of the 
temperature profiles that are obtained by hydrostatic equilibria in 
smooth potential wells; the not too peaked metallicity distribution; 
and finally the virial temperature of the gas (Tvir = 2.3 keV) that is 
not much lower than 3 keV (i.e., the declared limit of applicability 
of the Mazzotta et al.'s Tsi). A similar finding has been reported by 
Rasia et al. (2008). In fact, for models with strongly peaked metal- 
licity distributions, or with much lower mass (e.g., M = 10" M0 
and Tvir = 0.8 keV) we found that the two spectroscopic temper- 
atures are clearly different, and the Mazzotta et al. (2004) estimate 
would be higher than that of Vikhlinin (2006) (up to 20% in the 
explored To/T^ir range, for the quoted mass A/). 

Finally, it is useful to mention the relation between <T> and 
T500 (see end of Sect. 12. 2.4) . since the temperature profile is gener- 
ally recovered from observations out to radii smaller than r^ir, typi- 
cally out to rsoo with the most sensitive observations (e.g., Viklinin 
et al. 2006). The calculation of the T500/ <T> ratio for all our 
models confirmed the expectation that the hotter is the central re- 
gion with respect to the outer one, the higher is this ratio. In fact, 
we found for quasi-isothermal models Tsoo = (l-fl.5) <T>, 
as To /Tvir goes from 0.4 to 4; for quasi-polytropic models Tsoo = 
(1 1.4) <r>, as To /Tvir goes from 2.6 to 4. For TMB models, 
T500/ <T> varies respectively between 1.3 and 1.5, and between 
1.3 and 1.6, for a = 1 and a = 0, as rg/vc varies between 0.2 to 
2; for a = 2, Tsoo/ <T>= 1.3, the lowest value, a consequence 
of its cold central region. 



3.1.1 Quasi-isothermal and polytropic models 

We first discuss Fig. H For To>T^ir, <T>l and are both 
higher (up to a factor of ~ 2) than the mass-weighted temperature 
<T>, because they are dominated by the hotter central regions. For 



quasi-isothermal models with ro<Tvir, instead, the 3 temperatures 
almost coincide, because the gas is nearly isothermal (see Fig. 1). 
Overall, T^i and <T> l agree very well up to To ~ 2T^ir. Starting 
from ro>2rvir, <r> L is larger than T^i, a tendency that becomes 
stronger with increasing To /Tvir'. this is due to fact that Tsi is bi- 
ased towards the lower values of the range of temperatures [e.g. 
rai(3/4) weights each thermal component by p^T"''^^ instead of 
by p^A(r) ~ p^T^/^]. 

At high To the size of the discrepancy between <T> and < 
r> L or Tsi compares well with the analytical predictions based on 
the asymptotic profiles il4\ . i20\ . From these expressions, defining 
9 = <j)/<j)o, the limit values of <T> and Tsi (5), both in the quasi- 
isothermal and quasi-polytropic cases, are 



<T> = 



and 



T,,(5) = 



2k /(*- *t)'/^+*rf^x' 



(31) 



(32) 



which are independent of To. In particular, note that Tsi(2) corre- 
sponds to the case of pure bremsstrahlung emission, which is sim- 
ilar to the emission described by our adopted cooling function, at 
least for high temperature 

The analytical solution of the integrals l |31t and l |32t . for 
spherical 7 = 1 models with truncation at the virial radius, gives 
<r>~ l.llTvir, Tsi(2) ~ 2.29Tvir and T^i{3/4) ~ l.SOTvir. 
These values are close to those shown by <r> and <T>l at 
To ~ 4rvir (Fig.©; Tsi (3/4) instead is still far from its limit value, 



^ In the numerical computations we use the full expression for A(T, Z) in 
the 0.3-8 keV band, and this is not well represented by a simple power-law 



(see Appendix B). 



10 Ciotti and Pellegrini 



even though it has already started decreasing. As long as <r> l can 
be considered similar to rsi(2), then 7*81(3/4) is predicted to tend 
to~ 0.6 <T>L (which was verified with numerical models not 
presented here). The same calculations at the limit of high To for 
the 7 = potential give <T>~ l.lSTvir, rsi(2) ~ 2.27Tv,r 
and rsi(3/4) ~ l.SlTvir; these values are very similar to those for 
7 = 1, in agreement with the close location of the solid and dashed 
lines in Fig. [4] 



3.1.2 TMB models 

The TMB cases (Fig.O are more varied. The first result is that <T> 
(black lines) remains almost constant for different rg/rc, at a value 
nearly independent of a. Analytic integration for a = 2 shows that 
<T>= l.lSTvir, in agreement with the numerical result in Fig. 5. 
The second result is that again Tsi and <T> l overestimate <r>. 

Tsi (green lines) decreases steeply with increasing r^/rc for 
a = 0, because the density in the central hotter regions decreases 
by almost an order of magnitude, while it is increasing in the colder 
external region. This same trend is again present, though milder, 
for the models with a = 1. As discussed in Appendix A2, without 
a cut at low temperatures Tsi{S) diverges for the a — 2 models; 
therefore a cut at kT = 0.1 keV (which excludes the cold central 
regions where p oc r~^) has been adopted to produce Fig.|5] <T> l 
does not suffer from this problem, because of the low temperature 
cut in the cooling function. <T> l is higher than Tsi (for the same 
q) as in the quasi-isothermal and quasi-polytropic models, and for 
the same reason of being Tsi biased towards the lower values of the 
temperature range. 




r / r 



Figure 7. The mass weighted < Z > (black) and luminosity weighted 
<Z>£, (red) abundances calculated for TMB models within the 7 = 1 
potential, for a = (dotted lines), a = 1 (dashed lines) and a = 2 (solid 
lines), as a function of Tgjrc - For q = 2 the density and temperature pro- 
files are independent of Tg/rc, and so are <Z> and <Z>i^; for a = 0, 
only the values up to rg/rc = 1 are displayed, as for Fig.|5] 



3.1.3 Changing the dark mass amount and shape 

For all models (quasi-isothermal, polytropic and TMB) we investi- 
gated the effects of changing the total mass M and of flattening the 
dark mass distribution. We found that all the trends in Figs.|4]and|5] 
remain the same with a different total mass M and shape. 

For what is concerning the vahies of the average temperatures, 
<T>L /Tvir and Tsi/Tvir remain the same within few (~ 5) 
percent, for > 3 x W^^Mq. The only exception is Tsi/Tvir 
calculated according to Vikhlinin (2006) for the TMB models with 
Q = 2, that decreases by 13% going from M = 1.4 x 10^^ A/© 
to A/ = 3 X 10"A7q. <T> /Tvir and Tsi((5)/Tvir, instead, are 
independent of A4 for all models, as can be proved analytically: the 
curves in Figs. |4] and [5] depend (for all the other parameters fixed) 
only on To /Tvir or rg/rc, and on e and r;. 

For a fixed mass A/ (> 3 x we then changed the 

values of e and i] from zero to 0. 1 and 0.3, which produces a flat 
E7-like shape, and up to 0.5 and 0.5. The values of all the average 
temperatures, when rescaled for the different Tvir, remain the same 
within 5%. Again the largest variation is that of Tsi/Tvir calculated 
according to Vikhlinin (2006) for the TMB models with q = 2, 
that increases by 13% going from the spherical to the e = = 0.5 
shape (that corresponds to a very prolate ellipsoid). 



3.2 Abundance averages 

Figures [6] and [7] show the trend of <Z> and <Z>l as a function 
of To /Tvir (for quasi-isothermal and quasi-polytropic models) or 
fg/rc (for TMB models), for the same total mass A/ and potentials 
used in the previous figures, and for the abundance profile l l23b . 



with the parameter values specified below that equation. These fig- 
ures show that <Z> l (red lines) always overestimates <Z> (black 
lines), a result of the larger weight that central regions have in the 
calculation of <Z>l. For all models <Z> and <Z>l decrease 
with To/Tvir or rg/rc increasing; this is explained by the fixed 
metallicity profile coupled with gas density profiles that become 
flatter (Figs. [T] and |2j, so that the central regions where the abun- 
dance is highest become less and less important in the integrals of 
eqs. (1) and ( I27t (except for the TMB a = 2 models, where the 
density profile is independent of rg/rc). The steeper density pro- 
files are those of the isothermal models with To < 2Tvir (Fig. [T]!, 
therefore the decrease of <Z> and <Z>l at increasing To is more 
pronounced in this range of temperatures. 

In general, the details of the discrepancy between <Z> and 
<Z>L, and its trend with To/Tvir or Tgjr^, depend on how the 
density, temperature and abundance profiles differ from each other. 
The overestimate obtained by using luminosity weighted abun- 
dances is stronger for steeper gas density profiles. For example, 
the 7 = models have a slightly flatter density profile at the center 
than the 7 = 1 ones, so that the discrepancy between <Z> and 
<Z>L is in general slightly smaller for 7 = than for 7 = 1. 
A similar behavior is presented by TMB models. The profiles in 
this family may have a drop in temperature at the center, and con- 
sequently a steep increase in density, that is more pronounced for 
Q = 2 (Fig.[3j. Figure|7]reflects this fact, showing the largest over- 
estimate of < Z >L among TMB models for a — 2, while the 
smallest is that of a = models. 

In summary, considering all our models <Z>l / <Z> lies 
in a quite small range: 1.1< <Z>l / <Z> <1.6 for quasi- 
isothermal models, 1.25< <Z>l / <Z> <1.45 for quasi- 
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polytropic models, and 1.3< <Z>l / <Z> <1.7 for TMB 
models. 

3.2.1 Changing the dark mass amount and shape 

For all models presented in this paper it can be proved that <Z> 
is independent of M, while it depends on To/T^i^ or r^/rc, and 
e and As a consequence, the values of <Z> in Figs. [^^7] keep 
the same for all dark masses M when the metallicity distribution 
(23) is used with the specified parameters. <Z>i^ remains almost 
identical, for 71/ > 3 x IO^^A/q, with the largest variation for 
the highest To/T^ii and r^/rc of just 2%. These small variations 
are accounted for by the fact that the gas emissivity is not a pure 
power-law in temperature. 

As for the temperatures, we also investigated the effect of flat- 
tening of the mass distribution. <Z> and <Z> l become smaller 
when increasing the flattening with respect to the spherical case, 
which is explained by more and more gas mass being displaced 
at larger distances from the center, where the abundance is lower. 
When changing the values of e and from zero to (0.5, 0.5), for the 
range of To /Tvir of Fig.[6]and of rg/rc of Fig.l?] <Z> decreases by 
<16%, and <Z>l by <14%. Smaller variations are obtained for 
more reasonable flattenings, as for example of the order of ^ 5% 
for e and i] equal to 0. 1 , 0.3. 



4 SUMMARY AND CONCLUSIONS 

In this work we have compared the values of mass and luminosity 
weighted metallicity and temperature, for a large set of hydrostatic 
gas distributions, some of which resemble those typical of the in- 
tracluster and intragroup media. In addition, we also computed the 
temperature that would be derived from observed spectra (the so- 
called spectroscopic-like temperature) by using two recently pro- 
posed methods for its estimate. The results of this analysis are use- 
ful for distant groups/clusters, or in general for systems with a low 
number of observed counts, where only global average values can 
be recovered from observations. 

This study is based on a few steps. First, the potential well 
of triaxial dark matter halos, with different density slopes and ad- 
justable flattenings, was built analytically by means of homeoidal 
expansion, which gives a simple yet accurate analytical approxi- 
mation of the true potential. In the second step we showed how to 
construct hydrostatic analytical solutions for triaxial truncated den- 
sity distributions, and presented the equilibrium configurations in 
the quasi-isothermal and quasi-polytropic cases, and for a family 
of modified f3 models. In the third step we superimposed a metal- 
licity distribution derived from observations of the ICM/IGM. Fi- 
nally, the gas radiative properties were computed by using the cool- 
ing function A(r, Z) appropriate for our range of gas temperatures 
and a chosen sensitivity band of 0.3-8 keV. Mass and luminosity 
weighted temperature and abundances for the models were then 
obtained, thanks to the Projection Theorem. 

The main results can be summarized as follows. 

• The quasi-isothermal and polytropic models show gas density 
and temperature profiles similar to those observed for non-cool core 
clusters, while those of TMB models with a = 1 or q = 2 resem- 
ble cool-core clusters. In particular the temperature profiles of the 
latter TMB models, when rescaled to Tsoo , compare well in shape 
and normalization with observed profiles. In general, <T> /T500 
ranges between 1 and 1 .5, and it is 1.3 for TMB models with a — 2. 



• The luminosity- weighted temperature <T>l overestimates 
<T> up to a factor of ~ 2, and the discrepancies increase with in- 
creasing gas temperature (scaled by Tvir) for quasi-isothermal and 
polytropic models, or for increasing rg/rc for TMB models. For 
these latter models with a — 2 the overestimate is milder (a factor 
of ~ 1.3). 

• Tsi always provides a less serious overestimate of <T> than 
<T>L. The discrepancy bewteen T^i and <T> becomes smaller 
for increasing To /Tvir and for increasing rg/rc. 

• The exception to a general overestimate of < T > is that 
of "cold" (ro<1.2rvir) quasi-isothermal models, where the three 
temperatures <T>, <T>l and T^i are very similar. Also, <T>l 
and Tsi keep close up to To ~ 2Tvir, and depart for higher To /Tvir. 

• When changing the total dark mass M, the general behavior 
of <r>, <T>L and T^i described above remains the same. The 
values of <T>l /T^ir and Tsi/Tvir turn out to keep within 5% by 
changing M, for the range of masses typical of large groups and 
clusters (M > 3 x lO^M©). <r> /Tvir is instead independent 
of M. 

• In the explored range of triaxiality, flattening effects are not 
strong: the average temperatures normalized to Tvir remain the 
same within 5%. 

• The only exception to the small (<5%) variance with a change 
of shape or mass is given by the "cool-core" models (TMB models 
with a — 2): the increase of Tai/Tvir can be as large as 13% going 
from M = 3 X lO" Mq to M = 1.4 x IO^^Mq, or from spherical 
to e = 77 = 0.5 at fixed M. 

• The luminosity weighted < Z >l overestimates the mass 
weighted average abundance <Z>. For quasi-polytropic and quasi- 
isothermal models with To > 2.6rvir we found that 1.3< < 
Z >L I < Z > <1.5. This ratio extends over a larger range 
(1.1< <Z>L I <Z> <1. 6) for colder quasi-isothermal models 
(ro<2Tvir). TMB models show their smallest overestimate (~ 1.4) 
for a = 0, and the largest (~ 1.7) for a — 2, that repoduces the 
case of "cool core" ICM/IGM. 

• Similarly to what found for <T>, an M variation has no effect 
on <Z>, and a negligible effect on <Z>l. The effect of flattening 
is present, but it is not very important. For (e, r\) equal to (0.1, 0.3) 
and (0.5, 0.5), <Z> and <Z>l decrease by <few %, and by 
< 13% respectively. 

Thus, we have shown that when deprojection is not feasible or 
robust (as in the case of distant objects, significant deviations from 
spherical symmetry, etc.), the alternative approach of considering 
the global average values of temperature and abundance, obtained 
as surface integrals over the image, has the advantage over depro- 
jection of being independent of the shape of the system and of the 
relative orientation to the observer, but in presence of non-uniform 
metallicity and temperature distributions it must be calibrated by 
computing the appropriate correcting factors, as those determined 
in this paper. It would be interesting to apply both methods (depro- 
jection vs. surface average) to real systems with detailed observa- 
tions, and to compare the results. 
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APPENDIX A: ANALYTICAL RESULTS 
Al The Projection Theorem 

The numerical integrations performed in this paper are based on 
a very simple but far-reaching mathematical identity holding be- 
tween volume and surface integrals of projected properties of 
"transparent" systems of general shape (e.g., see Ciotti 2000). In 
the present context, let us consider a property P(x) (for example 
the ICM temperature or metallicity), associated with a field /^(x) 
(for example the ICM emissivity) acting as a "weight". Without 
loss of generality, we can suppose the line-of-sight to coincide with 
the 2:-axis, so (x, y) is the projection plane. The weighted projec- 
tion is naturally defined as 



T,{x, y) X Ppr(x, y)= j i/(x) P(x) dz 



(Al) 



where Ti{x,y) = u{x)dz. From the identity J E 

Ppr dxdy = J P(x) d^x, it follows that 



E X Ppr dx dy I i/(x) P{x)d''x 
<P>.= ^ — , (A2) 



E dx dy 



!/(x) d'^x 



where the integrals that appear in this definition extend on the 
whole image. In practice, the surface weighted integral of a pro- 
jected property coincides with the volume weigthed integral of this 
property; in turns, this proves that < P >v is independent of the 
viewing angle. 

A2 The density approach 

As well known, for a gas in hydrostatic equilibrium in a potential 
well a!)(x) 



(A3) 



and the density, temperature and pressure distributions are strati- 
fied on the equipotential surfaces, i.e., p = p{<j}), p = p(<?!>), and 
T = T[cj>), so that the gas is barotropic. In standard applications 
some relation p{p) is assigned, and eq. (A3) is solved. In the den- 
sity approach instead the functional form p{<j>) is assigned. In this 
case it is trivial to prove that the function 



H{<f>) 



p{t) dt 



(A4) 



satisfies the identity V// = -~p[(j>)V4>, and so integrating eq. (A3) 
over an arbitrary path it follows that 



(A5) 



where p, and are the pressure and potential values at some refer- 
ence position X,. For untruncated models in halos of finite mass the 
natural choice is to set p, = and (j)* = at x, = oo. For models 
as those considered in Sect. 2.2, with a finite value of the central 
potential (j>o and a truncated density p — pt, where pt = p{(f)t), it 
is natural to set p^ = pt — 0, so that eq. (A5) becomes 



p(x) = \<l,o\ - H(^t) - Pt X (* - *t)] 



(A6) 



where >I> = cp/(f>o and = p(*)d*. 

The simple idea behind the density approach is to choose a 
spherical gas density distribution p{r) with the desired radial pro- 
file, and a spherical potential <^(r). Due to the monotonicity of the 



potential (guaranteed by Gauss theorem), it is always possible (in 
principle) to eliminate the radius between the two distributions, 
thus obtaining p{<j}), and the function H. In the second step one 
deform the spherical potential (for example by using homeoidal 
expansion as in this paper, or with the complex shift method de- 
scribed in Ciotti & Giampieri 2007, or by functional substitution as 
in the Miyamoto-Nagai 1980 disk case) but assumes that the func- 
tion p{(l)) is the same as in the spherical case. As a consequence, the 
function H is still exact, while the deformed density profile is more 
and more similar to the spherical initial distribution for smaller and 
smaller potential deformations. 

In the family of models presented in Sect. 2.2.3, remarkably 
simple explicit cases are obtained for /3 = 2/3 and q = 0, 1, 2. In 
particular, after defining 9 as in equation (15), for a = we obtain 



Po 



6^* 



+ 



6(1 - b^) 



arctan- 



62 + 1 (62 + 1)2 
62 log[(l-*)2+62*2 



_6*_ 

L - * 



+ 



(62 + 1)2 



(A7) 



for Q = 1 



Hi^) ^^ 1- v/(l-»)2 + 62^2 



po 



62 + 1 



6l/» + ^(1 - *)2 + 62*2 t^(2 + b^) 
+ (62 + 1)3/2 X 



1 - (1 + 62)4' + Vl + 62,/(l-*)2 + 62*2 
log — — 7==^ (A8) 



1 + 



and finally 



(A9) 



Po 

for Q = 2. 

It can be easily proved, for example by series expansion and by 
order matching of the hydrostatic equation near the origin that the 
equilibrium temperature for r ^ converges to a positive value 
for a — 0, while for a > it vanishes as r" (0 < a < 1), as 
— rlogr (q = 1), and finally as r (a > 1). In particular, in the 
a — 2 models, the vanishing of temperature near the center is not 
strong enough to compensate the density (square) increase as r"**, 
and this leads to a divergent denominator in the definition of Tsi{5) 
for S < 5/2. 



APPENDIX B: THE COOLING FUNCTION 

The functions A(r, 0) (in units of 10"^'' erg cm^ s"^) and g{T) 
(dimensionless) appearing in eq. i25\ describe the cooling func- 
tion over the 0.3-8 keV band for plasmas in the temperature range 
kT — 0.1 — 16 keV and with the abundance ratios of Grevesse 
& Sauval (1998). In order to speed up the numerical integrations, 
in the code we use the following non-linear fits (the Fade approxi- 
mants) of the numerical values produced with APEC 

a + bT + cT^ + dT'-^ + eT* 



A(r,0) = 



f + gT + hT'^ + iT'^ + IT^ + mT^ ' 



and 



9iT) 



a + bT + cT'^ + dT^ + eT* + /T" + gT^ 



(Bl) 



(B2) 



h + iT + lT^ + mT^ + nT* + oT^ + pT^ ' 

In the expressions above T is expressed in keV, and the coefficients 
are given in Table B 1. 



14 Ciotti and Pellegrini 

Table Bl. Double-precision coefficients for the numerical cooling function 





A(T, 0) 


9{T) 


0, 


0.006 1 82527 1 25052084 


70.36553428793282 


b 


-0.19660367102929072 


-445.6996196648889 


c 


2.0971400804256 


1154.002954067754 


d 


-8.007730363728133 


-334.98255766589193 


e 


1 1 .575056847530487 


232.9431532001086 


f 


-0.005867677028323465 


185.87403547238424 


9 


0.1915050008268439 


-6.733225934810589 


h 


-0.8925516119292456 


-0.5099502967624234 


i 


1.2939305724619414 


39.52256435778057 


I 


0.4960145169989057 


-357.9367436666503 


m 


0.010460898901591362 


1343.6559752815178 


n 




-1877.1810281439414 


o 




1208.8015587434686 


P 




-27.997494448025837 



